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We evaluate Generalized Parton Distributions of the pion in two chiral quark models: the Spectral 
Quark Model and the Nambu-Jona-Lasinio model with a Pauli-Villars regularization. We proceed 
by the evaluation of double distributions through the use of a manifestly covariant calculation based 
on the a representation of propagators. As a result polynomiality is incorporated automatically 
and calculations become simple. In addition, positivity and normalization constraints, sum rules 
and soft pion theorems are fulfilled. We obtain explicit formulas, holding at the low-energy quark- 
model scale. The expressions exhibit no factorization in the i-dependence. The QCD evolution of 
those parton distributions is carried out to experimentally or lattice accessible scales. We argue for 
the need of evolution by comparing the Parton Distribution Function and the Parton Distribution 
Amplitude of the pion to the available experimental and lattice data, and confirm that the quark- 
model scale is low, about 320 MeV. 
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I. INTRODUCTION 

Generalized Parton Distributions (GPD's) encode de- 
tailed dynamical information on the internal structure of 
hadrons and have thus become in recent times a ma- 
jor theoretical and experimental endeavor (for exten- 
sive reviews see e.g. [3, [H, [3, 0, S IS 0, H| and refer- 
ences therein). They represent a natural interpolation 
between form factors and quark distribution functions. 
Actually, while form factors and distribution functions 
provide in a separate way the spatial and momentum 
quark distributions in a hadron respectively, GPD's pro- 
vide a simultaneous phase-space description of the quark 
hadron content as far as the position-momentum uncer- 
tainty relations allow Experimentally, GPD's show 
up in hard exclusive processes such as Deeply Virtual 
Compton Scattering (DVCS) or hard electroproduction 
of mesons. Factorization for hard exclusive electropro- 
duction of mesons in QCD was proved in Ref . [l(| ■ Ef- 
fects of the Regge exchanges to exclusive processes were 
investigated in Ref. [Til ] . 

In the present paper we are interested in GPD's for 
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the pion. Although there are little chances of mea- 
suring them directly in experiment, the pion GPD's 
are amenable to indirect experimental determination as 
well as studies both on transverse [l2j as well as Eu- 
clidean [l3j lattices (for a combination of experimental 
and lattice-based reconstruction in the proton case see, 
e.g., Ref. [14]). In addition, there are many theoreti- 
cal advantages for studying this quark- antiquark bound 
state. In the first place, in the chiral limit where the cur- 
rent quark masses vanish, the spontaneous breakdown 
of chiral symmetry of the QCD vacuum generates pions 
as the zero modes. Their properties are expected to be 
dominated by the broken chiral symmetry while confine- 
ment effects are expected not to be crucial. Actually 
finite mass corrections to GPD's of the pion have been 
treated in the standard [lj| and partially-quenched [l6| 
chiral perturbation theory, while the breakdown of the 
expansion for small x ~ m 2 /(47r/) 2 has been pointed 
out in Ref. 

[HI- 

Besides, compared to the nucleon there 
are no spin complications for the pion case, and thus the 
study reduces to two single scalar GPD function, one for 
each isospin combination. Finally, the pion provides a 
useful framework to learn on the interplay between the 
chiral symmetry and the light cone features, since we are 
studying the behavior of the would-be Goldstone boson 
in the infinite momentum frame. 

Despite the intrinsic complexity of the GPD's, there is 
a number of simple conditions which ultimately are con- 
sequences of the Poincare and electromagnetic gauge in- 
variance and provide a priori tests on the validity of the- 
oretical calculations. Proper support and polynomiality 
restrictions on the GPD moments \W are manifestations 
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of the Lorentz invariance. We note that polynomiality 
is not satisfied in light-front calculations. Double distri- 
butions do not suffer from this problem [l8| (using the 
double distributions is a way of projecting the Lorentz- 
violating term onto the right space) although they re- 
quire the so-called £>-terms [19( to comply with the most 
general polynomial allowed by the dimensional analysis. 
Normalization conditions and sum rules are a manifesta- 
tion of the gauge invariance, which at the quark level re- 
quires the correct implementation of the electromagnetic 
Ward-Takahashi identities. The positivity bound [20l.l2l| 
underlines the Hilbcrt-space quantum-mechanical prob- 
abilistic nature of pion light-cone wave functions, and 
may impose relevant constraints on admissible regular- 
izations based mainly on subtractions of ultraviolet di- 
vergences. Soft pion theorems based on PC AC relate 
GPD's to Parton Distribution Amplitudes (PDA's) [II]. 
On a theoretical level, the amazing aspect of GPD's is 
that the constraints that ought to be fulfilled a priori are 
so demanding and intricate that it is extremely difficult 
to provide ansatze fulfilling all of them simultaneously. 
This is why dynamical calculations going beyond reason- 
able but admittedly ad hoc parameterizations are partic- 
ularly interesting and instructive. On the other hand, 
dynamical models providing GPD's are also generating 
mutually consistent parton distribution functions, par- 
ton distribution amplitudes, and form factors. Although 
this may appear a rather trivial statement, it imposes 
demanding and tight constraints on details of the calcu- 
lation, and more specifically on the proper handling of 
ultraviolet divergences based on the correct implementa- 
tion of electromagnetic and chiral Ward-Takahashi iden- 
tities. Even in the case of a simple hadron such as the 
pion in the chiral limit, the above mentioned necessary 
conditions provide powerful limitations and in some cases 
clash with well established prejudices about the meaning 
and realization of a relativistic bound-state wave func- 
tions in quantum field theory [23j, l24l . l25j . 

In the present work we determine GPD's incorporat- 
ing all the desirable properties required by the sym- 
metries in two chiral quark models, the Nambu-Jona- 
Lasinio (NJL) (for reviews see e.g. H^,[23| and references 
therein) and the Spectral Quark model (SQM) [HH^, 
which essentially is a way of introducing regularization 
in such a way that the vector-meson dominance is built 
in. These field theoretical models incorporate the sponta- 
neous breaking of chiral symmetry. Chiral quark models 
use the large- N c expansion at leading order, where the 
observables are obtained with one-loop quark diagrams. 
The present calculation extends previous calculations of 
PDF's [HHlijl and PDA ' S El HI- Diagonal GPD in 
impact parameter space in these models were considered 
[321] . Our present GPD result reproduces consistently all 
these particular cases. Recently, the Transition Distribu- 
tion Amplitude (TDA) [H, [34| has also been evaluated 
in SQM [H. 

There has been a number of calculations of GPD's and 
related quantities of the pion within the framework in- 



corporating chiral symmetry. Early calculations of pion 
GPD's were done in an instanton-inspired model char- 
acterized by a momentum dependent mass function of 
a dipole form [3|| [37|. while PDF's were evaluated 
in the same model in [38l [39j . The crossing-related 
27rGPA was also evaluated in Ref. in the same in- 
stanton model disregarding the momentum dependence 
of the quark mass, a valid assumption in the limit of 
small instantons. In that limit, end point discontinu- 
ities arise. A full consideration of poles in the complex 
plane in a non-local version of the NJL model was de- 
scribed in Ref. ,4lj. Generally, the nonlocality of the 
quark mass function generated incorrect normalization, 
since as noted later, PC AC should be properly incor- 
porated [42| . an issue also emphasized more recently in 
Ref. [43d . A rather interesting feature of Ref. [12] is 
that end-point discontinuities reappear after PCAC is in- 
corporated, even for momentum-dependent quark mass 
functions, against the widely spread prejudice that they 
only arise for momentum-independent masses. The OPE 
and duality aspects of GPD's have been discussed in 
Ref. [44fl. Light- front calculations have been undertaken 
in Ref. [45l . |46| for point couplings with subsequent inser- 
tion of Gaussian wave functions, however the approach 
violates polynomiality. Power-law wave functions and 
GPD's of the pion proposed in Ref. [47| satisfy polyno- 
miality but violate positivity (see [481] V Studies paying 
particular attention to p olynomiality were first made by 
Tiburzi and Miller [48|, [4!| who proceeded via double 
distributions [501 ]. However, regularization in these works 
was done by introducing momentum dependent form fac- 
tors, which makes them difficult to reconcile with the 
gauge invariance. The model of Ref. [5l[ based on a pseu- 
doscalar pion-quark coupling does not incorporate chiral 
symmetry and does not fulfill the momentum sum rule. 
Nogucra, Theufil, and Vcnto carried out a calculation 
in the NJL model based on the light-front coordinates, 
where the fulfillment of polynomiality for non-vanishing 
momentum transfer is not apparent analytically [52j , and 
in fact numerical integration was required to establish 
this property. Our NJL results agree with that work, 
with the important methodological difference that the 
double distributions, where polynomiality is manifest, 
are used throughout. Moreover, our regularization is 
somewhat different than in the model of Ref. [521 ]. we 
also use the non-linear rather than linear realization of 
the chiral field. 

Despite the numerous model calculations of the GPD's, 
it remains to date unclear what is their significance or im- 
pact on the interpretation of actual experiments and/or 
lattice data. This is, perhaps, why most calculations of 
the GPD's based on dynamical quark models and going 
beyond just phenomenological parameterizations do not 
address this issue. However, while experimental or lattice 
results generate scale-dependent GPD's, embodying the 
well established logarithmic scaling violations in QCD, 
it is notorious that models generally produce scale inde- 
pendent functions. Thus, quark models represent those 
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distributions at a given low energy scale. It is noteworthy 
that scaling violations can only be computed in the twist 
expansion order by order. For instance, for the structure 
functions F(x, Q) with the Bjorken x and momentum Q 
one has for the quark model 



F(x,Q) = F (x) + ^- + 



while for QCD 

F(x,Q) = F (x,a(Q 2 )) + F2iX > am) 



Q 2 



(l.i) 



(1.2) 



where F n (x, a{Q 2 )) are low energy matrix elements with 
computable anomalous dimensions and depending log- 
arithmically on the scale through the running coupling 
constant 



a(Q 2 ) 



-in 



/? log(Q 2 /A2 ) 



0o = ^N c - 2 -N f . 



In this work we take 



Aqcd = 226 MeV 



(1.3) 
(1.4) 

(1.5) 



and N c = Nf = 3. The matching conditions are taking 
at a given scale Qo order by order in the twist expansion 



-Fn 0*01 Model = F n (x, a(Qo))|QCD- 



(1.6) 



A quite different issue is the operational definition of the 
low-energy reference scale Qn- Here we will use along the 
lines of previous works [231, H3, HH the momentum 
fraction carried by the valence quarks. It turns out that 
in order to describe the available pion phenomenology 
the initial scale Qo from the quark model must be very 
low, around 320 MeV. At such low scale the perturbative 
expansion parameter in the evolution equations is large, 
a(Qo)/2Tr — 0.34, which makes the evolution very fast 
for scales close to the initial value Qo- 

None of the previous chiral-quark-model studies of the 
genuine GPD's (off- forward non-diagonal) carried out the 
QCD evolution, starting from the initial condition pro- 
vided by the models. The evolution is a major element 
of this work. As already mentioned, it is also a crucial 
element if one wishes to compare the model prediction 
to the data from experiments or lattice simulations. At 
the moment these data are available only for the forward 
diagonal parton distribution function of the pion (PDF), 
or the PDA. 

We have taken an effort to separate formal aspects 
of the calculation from the model-dependent technical- 
ities of the regularization. That way we simply achieve 
the desired features on general grounds, such as the sum 
rules or the polynomiality conditions [l(. We stress this 
is achieved without a factorized form in the t variable. In 
one of the considered models (the Spectral Quark Model) 



the final results for the GPD's can be written in terms of 
rather simple but non-trivial analytic formulas, which al- 
lows for more insight into their properties. We also show 
that our GPD's satisfy the positivity bounds. In essence, 
all the known consistency conditions and constraints are 
indeed satisfied in our calculation. 

The outline of the paper is as follows: In Sect.[TTlwe list 
the definitions and properties of the GPD's. We intro- 
duce both the asymmetric and symmetric kinematics, as 
well as define the isospin projections. We give the generic 
quark-model expressions in terms of the two- and three- 
point functions. Sections IIIII and IIVI contain the results 
of the Spectral Quark Model and the NJL model, respec- 
tively. We discuss the general need for the QCD evolu- 
tion of chiral quark models in Section [Vl where we define 
the matching condition in the light of phenomenological 
analyses, as well as Euclidean and transverse lattice cal- 
culations. In particular, we show the evolved forward di- 
agonal PDF of the pion and confront it with the results of 
the E615 experiment at Fermilab 53]. This agreement is 
quite remarkable, but sets the quark model momentum 
scale to very low values, Qo — 320 MeV. Likewise we 
also discuss the evolved PDA as compared to the E791 
measurement |54| of the pion light-cone wave function 
and to the lattice data. We show how the QCD evo- 
lution leads to vanishing of the PDF at x = 1 and of 
PDA at x — 0, 1, which is the desired end-point behav- 
ior. The LO QCD evolution of our genuine GPD's, based 
on the standard ERBL-DGLAP equations, is carried out 
in Sect. IVI1 where the obtained quark- model initial con- 
ditions are evolved to higher momentum scales. Finally, 
in Section IVIII we draw our main conclusions. The Ap- 
pendices contain the technique of analyzing the GPD's 
through the use of the a representation for the propa- 
gators. This method, first introduced in the context of 
structure functions in Refs. [38l . |39| . allows for a man- 
ifestly covariant calculation and leads to simple formal 
expressions for the basic two- and three-point functions 
emerging in the analysis. The appearance of £>-terms is 
manifest and natural in this treatment, based solely on 
the Feynman diagrams in a conventional way. We also 
list explicitly the basic two- and three-point functions of 
the two considered model. We proceed via the double 
distributions, which leads to a simple proof of polynomi- 
ality 



II. DEFINITIONS AND QUARK-MODEL 
EXPRESSIONS 



A. Formalism for pion GPD 

The kinematics of the process and the assignment of 
momenta (in the asymmetric way) is displayed in Figs. [1] 
and [2 For the pions on the mass shell we have, adopting 
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the standard notation 



with the symmetry properties about the X — point, 



2 2 2 

P = m„, q 



-2p-q = t, 
0, p-n=l, q-n=-(. (2.1) 



Note the sign convention for t, positive in the physical 
region. 

The leading-twist off-forward [t ^ 0) non-diagonal 
(C 7^ 0) generalized parton distribution (GPD) of the 
pion is defined as 



(2.2) 



(n b (p + 9 )|^(0)7 • nT V(z)|7r a (p))| 2+=0 zi=Q , 

where < £ and the cc variable, — 1+C < a; < 1 , is defined 
in the asymmetric notation (c/. Fig. [I]), a and 6 are 
isospin indices for the pion, T is the isospin matrix equal 
1 for the isoscalar and T3 for the isovector cases, finally 
ip is the quark field and z is the light-cone coordinate. 
Explicitly, the two isospin projections are equal to 



dz 

47T 



(7r b (p + g)|^(0) 7 • n^(z)\7T a (p))\ z+=Q z±=0 . 



(2.3) 



(2.4) 



(7T b (p + 9)^(0)7 • r^(z) r 3 k») | z+=0 ^ =0 • 

One can form the combinations termed the quark and 
antiquark GPD's of the pion, 



(Hi=o(x,t,t)-Hi = i(x,t,t)). (2.5) 



From the general formulation it follows that H q (x,(,t) 
has the support x € [0,1], whereas Hq(x,Q,t) the sup- 
port x G [— 1 + C) C]- The range x € [0, C] is called the 
ERBL region, while x E + 0] and x € [C, 1] are the 
DGLAP regions, where the nomenclature refers to the 
QCD evolution, sec Sect. [VTl 

In the symmetric notation, somewhat more convenient 
in certain applications 1 , one introduces 



c 



x = 



x-C/2 



2-C' " l-C/2' 
where < £ < 1 and -1 < X < 1. Then 



H I=a ' 1 {X,^t)=H 



7=0,1 



e + x 2^ 



(2.6) 



(2.7) 



1 In this paper we switch back and forth between the two conven- 
tions, since explicit expressions are shorter in the asymmetric 
notation, while some formal features are simpler to state in the 
symmetric notation. 



H^iX^t) = -H I=0 (-X,£,t), 
H I=1 (X,^t) = H I=1 (-X,£,t). (2.8) 



The following sum rules hold: 



JdXH I = 1 (X,Z,t) = 2F v (t), (2.9) 



dx x h i=0 (x, t) = e 2 (t) - eoi(t), (2.10) 



where Fy(t) is the electromagnetic form factor, while 
6*i (t) and 62 (t) are the gravitational form factors of the 
pion (see Appendix [Dj which satisfy the low energy theo- 
rem 6>i(0) = 6» 2 (0) in the chiral limit [55]. Sum r ule g 9j) 
expresses the electric charge conservation, while f|2 . 10|) is 
responsible for the momentum sum rule in deep inelastic 
scattering. Finally, for X > 

H I =°> 1 (X,0,0) = q(X), 

relating the distributions to the the pion's forward diag- 
onal parton distribution function (PDF), q(X). 
The polynomiality conditions [J state that 

f dX X 2j H'^iX, U)=J2 A i } 

i=0 

f l J'+l 

dX X 2 i +1 H I=0 (X,Z,t) = ^B t 0) (t)e 21 , (2.11) 

i=0 



where A [ f (t) and b\ 5) {t) are the coefficient functions 
(form factors) depending on j and i. The polynomiality 
conditions follow from basic field-theoretic assumptions 
such as the Lorentz invariance, time reversal, and her- 
miticity, hence are automatically satisfied in approaches 
that obey these requirements. Conditions (|2. 1 1[) supply 
important tests of consistency. In our approach the poly- 
nomiality will be demonstrated straightforwardly in an 
analytic way through the use of double distributions, see 
Appendix IA"1 

Another constraint for the GPD's, the positivity 
bound [2l| , is derived with the help of the Schwartz in- 
equality and the momentum representation of the pion 
light-cone wave functions. In the simplest form the con- 
straint states that (for t < 0) 



\H q (X,£,t)\ < ^/q{x in ) q (x out ), £<X< 1.(2.12) 

where x in = (z + £)/(l +£), x out = (x - £)/(l - £). 

The off-forward (A ± ^ 0) diagonal (£ = 0) GPD of 
the pion (we take 7r + ) can be written as 



H(x,£ = 0,-Al) = / d 2 be 



lAi'b 



q(x,b). (2.13) 
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We use here 

q(x,b) = J ^e^ +z ~ (2.14) 

x (w + (p')\q(0, -£:,b)7 + g(0, ^,b)|7r+(p)), 

where a; is the Bjorken x, A± = p' — p lies in the trans- 
verse plane, and b is an impact parameter. The model- 
independent relation found in Ref. (56| reads in the pion 
case 

j\xq(x,b) = J ^e^Fvi-ql) (2.15) 

By crossing, the process related to the Deeply Virtual 
Compton Scattering (DVCS) off the pion, i.e., two pion 
production in 7*7 collisions, can be measured at low in- 
variant masses [57|. The relevant matrix element reads 

<5> ab (u,(,W 2 ) = [ — e ,lp+r x (2.16) 
J 4tt 

(K a { Pl )K\p2)\mi ■ nTiP(z)\0)\ z+=Q z±=Q , 

where W 2 = {P1+P2) 2 , C = Pi-n/P-n and u — {p\—p2) 2 - 
By comparing, we have 

<S> ab {uX,W 2 ) =H ab (x,C,t) (2.17) 

One has the soft pion theorem [22] , 

$ /=1 (u, 1, 0) = H I=1 (2u - 1, 1, 0) = ((>(u), (2.18) 

where 4>(u) represents the Pion Distribution Amplitude 
(PDA) defined as 

(7r a (p)|^(z) 7 ^ 75 irV(0)|0)| z+=0i2i=0 

= ifp»6 ab f dxe txp - z <t>{x). (2.19) 

Note that result (|2.18|) is based on PC AC and hence is a 
consequence of the chiral symmetry. One of the reasons 
to prefer GPD's rather than 27rPDA is the absence of 
final state interactions, which are suppressed in the large 
N c limit 2 . 




p k-p P+q 




p p+q 



FIG. 1: The direct (a) and crossed (b) Feynman diagrams for 
the quark-model evaluation of the GPD of the pion. 



B. Formal results for chiral quark models 

The reduction formulas applied to the definition (|2.4|) 
result in the amputated three-point Green function with 
the constrained quark momentum integration, k + — xp + . 
Large-A^c treatment leads to one-quark-loop diagrams, 
with massive quarks due to spontaneous chiral symmetry 
breaking. In nonlinear chiral quark models the quark- 
pion interaction is described by the term — r ijjujU 5 ip in 
the effective action, where the pion field matrix is 

U 5 = exp(j 75 r • 0//), (2.20) 

where/ denotes the pion decay constant. The resulting 
Feynman rules and the definition (|2.2[) lead to the Feyn- 
man diagrams of Figs. [T] and f2p. The presence of the 
contact term is crucial for the preservation of the chi- 
ral symmetry [19]. The evaluation of the diagrams is 
straightforward, giving the following result for the isos- 
inglet and isovector parts: 

H I=0 (xX,t) = H a (x,(,t)+H b (x,(,t)+H c (x,(,t), 
Hi =1 {x,C,t) = W«(as,C s t) -H b (x,C,t). (2.21) 



2 The simplest example illustrating this feature is provided by the 
pion electromagnetic form factor. The radius reads 

4N C 



v 



1 



■ log 



Ml 



the first contribution stemming from the quark loop and the 
second contribution an estimate from pion loops |58l . 



The explicit contributions of the subsequent diagrams 



3 The similar calculation of Ref. uses the linear realization of 
the chiral symmetry, with the a field present. 
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responsible for the D-term. 



to the GPD's of Eq. (f2T2T|) are 

iN r uj 2 



Arr 2 f 2 



d kS(k ■ n — x) x 



l -k 2 -C, (uj 2 -k 2 + k -p)+x (uj 2 -^-k 2 + 2k -p)-k ■ q 
D k D k , n D, 



k-Uk+q^k-p 



t \ (Fk5{k -n-x)x 



-LU 2 + k 2 + (k ■ p+x {u 2 - 5-fc 2 -2k ■ p-2k ■ q) +k ■ q 

DkDk+qDk+p+q 

n o(x, (, t) = ^! f d 4 kS ( k ■ n - x) C . (2.22) 



4tt 2 / 2 



DkDi, 



k^k+q 



The denominator of the proparator of quark of mass to 
and momentum / is denoted as 



Di = l 2 -u 2 + iO. 



(2.23) 



The powers of the momentum k in the numerators may 
be eliminated with the following reduction formulas: 



k 2 = D k +u 2 , 

1 

2 ( 



k-q = -(Dk +q -t- D k ), 



k-p = ~(D k _ p - m 2 - D k ), 

k-p = \(Dk +P+q -ml- D k +q + t). (2.24) 
Then the GPD's become 

n w (x, C, t) = ~ l ^f 2 I d 4 k5(k ■n-x)x (2.25) 



1 



8tt 2 / 2 

i-C 



DkDk-p Dk+qDk-p Dk+qDk+p+ q UklJk+p+q 

(C - 2a;)™ 2 + t(x - 1) (C - 2a;)™ 2 + t{x - ( + 1) 



DkD k 



DkDk+qDk- 



T 



DkDk+qDk+p+q 



with the upper (lower) signs corresponding to the case 
of / = (I = 1). Note that the piece with l/{D k D k + q ) 
cancels out due to the presence of contact diagram (c). 



FIG. 3: The contour C for evaluation of observables in the 
meson dominance variant of SQM. Mv denotes the p-meson 
mass. The cross and hatched regions indicate the position of 
the pole and cuts of the spectral function Eq. 



The contribution of the diagram (c) , having the support 
for x £ [0, C], is the D-term 

From the above form it is clear that we need to consider 
two generic types of two- and three-point integrals: 



I(x, I- n, I' ■ n, {I - I') 2 ) = 

J(x, I ■ n,l' ■ n, I 2 , 1' 2 , 1 ■ I') 
iN r uj 2 



-iN c cj 2 
4tt 2 / 2 



, 4 S(k ■ n - x) 



d"k 



At: 2 ! 2 



d , k S(k.n-x) 
D k Dk-iD k _i> 



D k -iD k _i> 



(2.26) 



These are analyzed in detail Appendices IA II and IA 21 
The two-point function / is logarithmically divergent, 
hence the analysis needs regularization. This is where 
different quark models depart from one another. We 
may separate the issues of regularization from formal ex- 
pressions, which is convenient for theoretical aspects and 
the demonstration of the consistency conditions. Writ- 
ten in terms of the basic two- and three-point functions 
Eq. (|2T23| become 

H I= o,i(x, C, t) = X - [I(x, 0, 1, ml) + (1 - QI(x, C, 1, ml) 
Tl(x, -1 + C, C, m 2 ) T (1 - Ql(x, -I + (, 0, ml) 
-[(C - 2a;)m 2 + t(x - l)}J[x, (, 1, t, ml , -~] (2.27) 

±[(C - 2x)m 2 + t(x -( + 1)] J[C - x, (, 1, t, ml, -I] 



This equation may be considered as the generic non- 
linear local chiral quark-model result for the isospin- 
projected GPD's of the pion. Model details, such as reg- 
ularization, affect the specific form of the two- and three- 
point functions, but leave the structure of Eq. (|2.27p un- 
changed. The nontrivial features of the regularization 
will utterly be responsible for the fulfillment of the gen- 
eral properties of GPD's described in Section Hi Al 
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III. RESULTS OF THE SPECTRAL QUARK 
MODEL 



Now we come to the evaluation of GPD in specific mod- 
els. From now on we work for simplicity in the chiral 
limit, 



My = m p = 770 MeV. The contour C for the integration 
in (|3.2p is shown in Fig. [3] Despite the rather unusual 
appearance of the spectral function, the model leads to 
conventional phenomenology [2^.[59|. Importantly, it im- 
plements the vector-meson dominance, yielding the pion 
electromagnetic form factor of the monopole form 



= 0. 



(3.1) 



The first model we consider is the Spectral Quark Model 
(SQM) of Ref . , where all the necessary details of the 
model can be found. The one-quark-loop action of this 
model is 



v (t) ~ Ml 



My 



For the gravitational form factor we find 



(3.4) 



iN c / dwp(w)Tr log (i0 - ujU 5 ) 



(3.2) 



where p(ui) is the quark generalized spectral function, and 
U 5 is given in Eq. (|2.20p . In the calculations of this paper 
we only need the vector part of the spectral function, 
which in the meson- dominance SQM [291 ] has the form 



pv{uj) 



1 1 1 

2^iu (l-4w 2 /Af 2 ) 5 /2 



(3.3) 



exhibiting the pole at the origin and cuts starting at 
±My/2, where My is the mass of the vector meson, 



My 



_ pSQM 

= fa 



(*)• 

(3.5) 



Both the electromagnetic and gravitational form factors 
for SQM are plotted in Fig. 2] with solid lines. 

With the help of Eq. (|B1IB8|) it is straightforward to 
obtain the formulas for the GPD's in SQM. The ex- 
pressions are simple in the chiral limit, and shortest in 
the asymmetric notation. For the quark and antiquark 
GPD's we obtain 



2H q (x, C, t) = 0((1 - x)x) + 0((1 - x)(x - 0) + t(l - x) 



2(x - 1) (t(x - l) 2 + 3(C - 1)M 2 ) 6»(1 - x)9(x - C) 
(t(x - l) 2 + (C - l)M v f 



(x - 1) (t(x - l) 2 + 3(C - 1)M 2 ) (x(( - 2) + C) (3(C - l)C 2 M v + t ((C 2 + 8C - 8) x 2 + 2(4 - 50O + C 2 )) 



{t{x-iy + {Q-l)M v f 
x 9(x)6((-x)}, 



V 



(c 2 + ^p 1 ) 3/2 (^-i) 2 + (c-i)M 2 f 



n q (x,c,t)=n q (c-x,c,t). 



(3.6) 



From these, the isospin combinations are trivial to 
get. The formulas satisfy the consistency relations 
(|2.8l2.9l2.10p . In particular, upon passing to the sym- 
metric notation and using the above formulas we verify 



ERBL and DGLAP regions: 



f dXH I=1 {X^.t) = 2Fy- QM {t), 

J" dX X H I=0 (X, e, t) = (1 - eWl QM (t). (3.7) J, dx H i {x > C t} 



The isovector norm is decomposed as follows between the 



2-C M v C(M v +t(l-Q) 
2 Ml -t(2-()(M v -t(l-C))' 



dxH q (x, (, t) — 



2-C Ml 



2(1 - C) (M v - t) 



2 M v -<(2-C)(M 2 -t(l-C))' 



(3.8) 
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Fv(t) 




-t [GeV^ 



-t [GeV^] 



FIG. 4: (Color online) The pion electromagnetic, Fv(t), (top) 
and gravitational, 8\(t) = 62(f) = Fs(t), (bottom) form fac- 
tors in SQM (solid line, Eq. Q3.4l3.5p ) and in NJL model 
(dashed line, Eq. (jO)l). 



Some special values of the GPD's in SQM are: 

W I=1 (l > C,t)=W /=0 (l,e,t) = l, (3.9) 



M 2 {M 2 +t(l-Q) 
(M*-t(l-0? 



lim H I=0 (x,C,t) 



M*(M*+t(\-Q) 



«-c+ (M 2-t(i-c)) 2 ' 

T n , s t (3Mf - i(l - 0) (1-0 
lim W /=0 (a;,C,t) = — , • 



The values at x = — 1 + ( and a; = follow from the 
symmetry relations H I=10 (x, (, t) = ±H I=1 ' Q (C,-x, Q, t). 
We note that the discontinuities at the end points x — 
±1 and for the 1 = part at x = and x = C are a 
typical feature of quark-model calculations. The QCD 
evolution immediately washes out these discontinuities, 
see Sect. [VI] The derivative of H I=1 (x, £, t) with respect 
to x is continuous at the point x = where 



2M^ (3Af2 + t(l - 0) 



(M^-t(l-C)r 



(3.10) 



On general grounds, Eqs. (13. 6|) also satisfy the poly- 
nomiality conditions (|2.1ip . which can be seen from the 



derivation through the double distributions shown in Ap- 
pendix El Also note that the obtained formulas are cer- 
tainly not of the form where the t-dependence is factor- 
ized, i.e. 



H 



1=0,1 



(X,t,t)?F(t)G(X,Z). 



(3-11) 



For the case of t = the formulas (|3.6[) simplify to the 
well-know [l^, [52j step-function results 

n I=0 (x, C, 0) - 6[(1 -x)(x- C)] - 0[-ic(a: + 1 - C)], 
n I=1 (x,C,0)=9[(l-x)(x + l-Q]. (3.12) 

Another simple case is for £ = and any value of i, 

M£ (M*+t(x-l)*) 



n q (x,o,t) 



(M*-t( X -m 2 



(3.13) 



which agrees with the result reported in [321 ]. The cor- 
responding impact parameter representation obtained 
there is given by the formula 4 



q(x,b) 



M 2 



2tt(1 - x) 2 
bM v 



(3.14) 



-Ki 



-Kq 



fbM v 



1 — x ± \1 — x J "\l-x 
From this expression one obtains 

M 2 K (bM v ) 



dx q(x, b) 



2tt 



(3.15) 



This complies to the model-independent relation (|2.15|) 
when the vector-dominance form factor (|3.4p is used, 
since, explicitly, 



d 2 q 



iq-b 



2tt M 2 



K Q (bM v ). 



(3.16) 



The case of Eq. <\3.6\i for £ = 0.5 and several values 
of t is shown in Fig. Results for other values of C are 
qualitatively similar. Figure [5] shows the isospin and 1 
combinations, Hi = o,i- We note that the 1=1 GPD and 
its first derivative with respect to x is continuous at x = 
and x = £, while the 1 = combination is discontinuous 
at these points. 

At t = we have the above-mentioned step-function 
results 

H I=1 = 6(1 - X 2 ), (3.17) 
H I=0 = 6((1 - X)(X - 0) - 0((l + X)(-X - 0). 

As —t increases, the strength moves to the vicinity of the 
X = ±1 points. The limit oi-t(l-x) -> 00 in Eqs. ([3~6]) 



Note an overall sign missing in Ref. [32 
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H (x,r,t 



1 . 5 



0.5 



-0.5 




FIG. 5: (Color online) The SQM results for the quark GPD 
of the pion, Tt q of Eq. (|2.5[1 . plotted as a function of x for 
C = 1/2 and t = 0.2, 0, -0.2, -1, -10, -100 GeV 2 , from top 
to bottom (at x — 0.9). Asymmetric notation. 



yield the asymptotic forms 



W J=1 > (l,C,t) = l, 



(3.18) 



H I=1 (x,C,t) 



M^(l-C) 
t{l-x) 2 ' 

M 2 (W) 



€ [0, 1) 



i(l — .t)^ 

M 2 

W /=o fr , rt „ M v(2^-C) (C 2 -3C + 2) 
n W"*'- 2(x-l) 2 (x-C + l) 2 t 



, ice(o,0- 



In the DGLAP region the absolute value of the / = 0, 1 
functions are bounded by unity. Note that at large —t 
the GPD's continue to be equal to 1 at x = 1, however 
very quickly drop to zero in the DGLAP region. In the 
ERBL region the 1=1 part drops, while the 1 = part 
tends to —1 as x — > f - , and drops to elsewhere. 

Since in our quark model the parton distribution func- 
tion is unity, q(x) = 1, the positivity bound (|2.12|) states 
that (for t < 0) 



\H^ M (X,^t)\<l, £<X<1. 



(3.19) 



It is a priori not obvious that the bound should hold in 
chiral quark models where finiteness of observables re- 
sults from regularization involving subtractions. Never- 
theless, we have checked with Eqs. (|3.6p that condition 
p,19p is actually satisfied in the DGLAP region for all 
values of £ and all negative t. This is also manifest in 
Fig. O as well as in Eq. (|3.18|) . The bound is saturated 
at the end points X = ±1. Thus the positivity bound is 
satisfied in SQM. 




H I = (X, f,t) 




-Q-.-5--- 



FIG. 6: (Color online) Same as Fig. [5] for Hi = i (top) and 
Hi = o (bottom) in the symmetric notation, £ = 1/3. 



IV. RESULTS OF THE 
NAMBU-JONA-LASINIO MODEL 



We use the non-linear NJL model with Pauli-Villars 
(PV) regularization in the twice-subtracted version of 
Ref. (27j . The prescription for regularizing an observable 
O in this model is 



rcg = 0(O)-0(A 2 ) + A : 



dO{A 2 ) 
dA 2 : 



(4.1) 



where A are the PV regulator. Note that Eq. (|4.ip is 
different from the prescription used in [52| , where a vari- 
ant of the PV regularization with two distinct cut-offs is 
applied. We also use the non-linear rather than linear 
realization of the chiral field. In what follows we take 
M = 280 MeV for the quark mass and A = 871 MeV, 
which yields / = 93.3 MeV [27| according the the formula 



3M 2 

47T 2 



(log(A 2 + M 2 )) i 



(4.2) 



The pion electromagnetic form factor in the NJL model 
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Hi = i(X,f,t) 
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FIG. 7: (Color online) Same as Fig. for for the NJL model 
with the PV regularization, M = 280 MeV, A = 871 MeV. 

is 

F^(t) = l + ^Jx (4.3) 

/2V 4(Af»+A»)-t 1og f VW^Z!^ ) \ 
v v ; 5 \ V 4 ( M2 + A2 )-*+V=* / 

V / rcg 

The property Xva\ t ^>-ooFm\Sfy — follows from 
Eq. (|4.2p . The isovector form factors arising in both 
considered models are compared in Fig. 2] The formula 
for the gravitational form factor in the NJL is lengthy, 
hence we only give the numerical results in Fig. |4j In this 
model also the two gravitational form factors are equal 
to each other, 0f JL (t) = 9^ JL {t). We note that although 
the form factors in both models are qualitatively simi- 
lar, they are quantitatively somewhat different, which is 
partly due to the choice of parameters in the NJL model, 
as well as follows from different analytic structure of the 
corresponding formulas, in particular at large values of 
-t. 

The application of the formulas derived in Appendix [Cl 
leads to expressions similar to those of [52j • The analogs 
of Eq. (|3.9I3.18[) in the NJL model are more complicated, 



hence we do not give them here. More details may be 
found in Ref. [52j . For the special case of t = Eq. (|3.12|) 
holds. As in SQM, the conditions l|2.9|2.10|2.8|) are sat- 
isfied in the considered NJL model with the PV regular- 
ization. 

The numerical results for the NJL model are displayed 
in Fig. [7j When comparing Figs. [6] and [7j we note a strik- 
ing similarity between the two considered quark models. 
The slight differences stem mainly from different form 
factors in the two considered models, cf. Fig.|4j Our re- 
sults are also qualitatively similar to the case of the chiral 
limit in Fig. 6 of Ref. [52^ . As pointed out in AppendixPCl 
polynomiality follows from the derivation proceeding via 
the double distributions. We have checked that the simi- 
larly to SQM, the positivity bound (|2.12[) is also satisfied 
in the NJL model at any £ and all negative values of t. 



V. QCD EVOLUTION OF QUARK MODELS 
A. The need for evolution 

A key question, not only for our model but for any non- 
perturbative calculation, is what is the scale at which our 
model result for the GPD's holds. Ultimately, this boils 
down to the issue on how the model predictions for the 
GPD's might be confronted to experimental or lattice 
data. In QCD, the GPD's are scale dependent, while in 
models they correspond to functions defined at a given 
scale. This is so because low energy models hold at a 
scale above which scaling should set in. Perturbative 
QCD and the corresponding scaling violations bring in 
the issue of evolution equations for GPD's which will be 
treated in detail in Sec. EH In this section we discuss and 
update the procedure already used in previous works (23l . 
123, Mill for the evolution of PDF and PDA and extract 
its consequences as compared to available experimental 
data or lattice results. 

From the point of view of perturbative QCD where 
both quarks and gluons contribute as explicit degrees of 
freedom, the role of the low-energy chiral quark models is 
to provide initial conditions for the QCD evolution equa- 
tions order by order in the twist expansion. Clearly, chi- 
ral quark models contain non-perturbative QCD features, 
particularly the spontaneous chiral symmetry breaking. 
On the other hand, chiral quark models do not contain 
the QCD degrees of freedom, i.e., the current quarks and 
explicit gluons. So one expects typical high-energy per- 
turbative QCD features, such as radiative corrections, to 
be absent in the model. This is precisely the pattern of 
logarithmic scaling violations which the models lack but 
which have traditionally been computed in the perturba- 
tion theory in QCD. 

The procedure applied in this paper takes the quark- 
model distributions at some low quark-model scale Qo 
and evolves them to higher scales, where (for some ob- 
servables) the experimental or lattice data are available. 
In the following we use the leading-order ERBL-DGLAP 
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evolution equations with three flavors. This strategy re- 
flects the present state of the art, which can be validated 
by comparing our predictions both to experiment as well 
as available lattice data. It should be noted, however, 
that to date there is no rigorous relation between the 
QCD quarks and the constituent quarks of the chiral 
models, and a more fundamental description of the tran- 
sition from the hard to the soft would be very helpful. 



B. Momentum fraction and the matching condition 

For definiteness, we consider 7r + , and denote q{x) and 
q(x) the single-flavor distributions of quarks and anti- 
quarks. The valence (or non-singlet) quark distribution 
is 

V = U n — Uir + d-K — (In, (5.1) 

while the non-singlet quark distribution is 

S = Ujt + u n + d-x + Jjt + + s v . (5.2) 
The sea quark distribution is defined as 

s = S — V = 2(uk + d„) + s n + s„. (5.3) 

Isospin and crossing symmetries implies the property 

u n +(x) = d„+(l — x). (5.4) 

The energy-momentum tensor Q a/3 is a conserved 
quantity in any relativistic theory and hence renormaliza- 
tion invariant, due to the Poincare invariance. Its diago- 
nal matrix element between the pion state of momentum 
p is 



?r(p) \@ ap \T:{p)) = 2p a p 



- Or/*™/ 3 



(5.5) 



For the QCD Lagrangian the energy-momentum tensor 
can be separated into several contributions in a gauge- 
invariant but scale- and hence scheme-dependent man- 
ner (6(|. Although we will be considering the LO evo- 
lution, for our purposes we may have the standard MS 
scheme in mind and write 



(5.6) 



where Qg 13 , O s , and O"' 3 are the gluon, sea-quark, and 
valence-quark contributions, respectively. They are equal 
to 

H6f = 2pV<* >„(/*) 
= 2pV<*>.(/i) 

(5.7) 

where (x) s (/x), (x) s (p) and (x} v (fi), are the gluon, sea 
quark, and valence quark momentum fractions of the pion 



at the scale respectively. In deep inelastic scatter- 
ing (DIS) it can be shown [6l| that if q(x,fi), q(x,[i) 
and G(x,/i) represent the probability density of finding a 
quark, antiquark, and gluon, respectively, with the mo- 
mentum fraction x at the scale fj, (typically, we identify 
fi 2 with Q 2 in DIS), then 



dxx G(x, n) 



(5.8) 



(x) s (p) = / dxxs(x,n) 
Jo 

(x) v ([i) = / dxx V(x,fx). (5.10) 



(5.9) 



where, due to the crossing symmetry (|5.4[) for a sin- 
gle flavor one has (x) q — (x) u = (x)j = (x) v /2. The 
scale-dependent momentum fractions fulfill the momen- 
tum sum rule 



,00 



.(/*) 



;(M) 



1. 



(5.11) 



which is a consequence of the energy-momentum tensor 
conservation. In perturbation theory due to radiative 
corrections (x) g (p) and {x) s (fx) decrease as the scale /x 
goes down. On the contrary, the valence contribution to 
the energy momentum tensor evolves as 



(x)v(Q) 



{x}v(Qo) \a(Qo 



a(Q) 



7i 0) /(2/3o) 



(5.12) 



where j[ 0) /{2f3 ) = 32/ 81 for 7V F = N c = 3. Down- 
ward LO QCD evolution would yield that for some given 
reference scale, fi Q = Q Ql 

(x }„ (Q ) = 1 (x) s (Q ) + (a; ) g (Qo) = 0. (5.13) 

The scale Qo defined with the above condition is called 
the quark model scale for obvious reasons, as only va- 
lence quarks contribute. This may represent the match- 
ing condition between QCD and the chiral quark models, 
schematically written as Eq. (|1.6p . 

There exists a wealth of information on the momen- 
tum fraction carried by valence quarks in the pion at 
scales Q ~ 2 GeV, coming from several sources. Phe- 
nomenological analyses require these high scales to ne- 
glect higher twist corrections. The Durham group [62l |. 
based mainly on the E615 Drell-Yan data [53| and the 
model assumption that sea quarks carry 10 — 20% of the 
momentum fraction, determines (x) q = 0.235(10) at the 
scale Q — 2 GeV. The analysis of Ref. [63j, based on 
the assumption that the momentum fraction carried by 
valence quarks in the pion coindides with that of the nu- 
cleon, yields (x) q = 0.2 at Q = 2 GeV. 

Other determinations, comprising lattice calculations, 
may access directly the leading-twist contribution in 
a non-perturbative manner. However, the transition 
from the intrinsically non-perturbative lattice regular- 
ization to the perturbative MS regularization scheme re- 
quires high-energy matching scales. Early Euclidean lat- 
tice simulations provided (x) q = 0.32(5) at the scale 



12 



Q 2 w 4.84 ± 2.2 GeV 2 [H]. More recently, lattice cal- 
culations linearly extrapolated to the chiral limit [65j . 
yielding (x) q = 0.28(1) at the scale Q 2 « 5.8 GeV 2 , a 
somewhat larger value than suggested by phenomenol- 
ogy [62j, |63| and expected from the quenched approxima- 
tion. Still, in the quenched approximation, in Ref. [66| 
(x) q = 0.243(21) at Q = 2 GeV for light pions, which 
is in a closer agreement to the Durham [62| than to the 
Dortmund [63| results. This value squares with the gluon 
content of the pion (x) g = 0.37 ± 0.08 sta t ± 0.12 sys at a 
similar scale but with ~ 900 MeV, as extracted re- 
cently in Ref. [671 ]. 

Finally, there exist transverse lattice calculations, 
where full x-dependent parton properties can be deter- 
mined non-perturbatively at low scales [12]. The calcu- 
lation of Ref. il gives (x) q = 0.43(1) at Q 2 ~ 1 GeV 2 , 
whereas Ref. [63] provides, still at very low scales Q 2 ~ 
0.4 GeV 2 , the value (x) q « 0.38. 

For definiteness we adopt the values used in previous 
works (H H3, |U S3, namely that at Q 2 = 4GeV 2 the 
valence quar ks carry 47% of the total momentum in the 



(5.14) 



xq(x) 



pion 62, 66], e.g., for 7r H 



(x) v = 0.47(2), 
at Q 2 — 4 GeV 2 . At LO the scale turns out to be 
Qo = 313t 2 >eV, 



(5.15) 



where the value of Aqcd is provided in Eq. (|1.5|l and 
the error reflects the uncertainty in Eq. (|5.14|) . At such 
low scale the perturbative expansion parameter in the 
evolution equations is large, a(Q 2 )/(2Tr) = 0.34, which 
makes the evolution very fast for the scales close to the 
initial value. We return to this issue below. 



Evolution of PDF 



In Fig. [8] we display the forward diagonal PDF of the 
pion, in particular the quantity xq(x) = xH q =1 (x, 0, 0) 
obtained with the LO QCD evolution up to Q = 4 GeV 
from the quark-model initial condition, q(x, Qo) = 1 [23| . 
At leading order the standard DGLAP evolution holds, 
which for the Mellin moments reads 



dxx n qM)= j_(?m_ 

o n+1 \a(Qo) 



7, ( ° ) /(2/9o) 



(5.16) 



where the anomalous dimensions for the vector vertex are 
given by 
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FIG. 8: (Color online) The quark model prediction for the 
valence parton distribution (PDF) of the pion for a single 
quark (either u or d for 7r + ) evolved to the scale of Q = 4 GeV 
(band). The width of the band indicates the uncertainty in 
the initial scale Qo, Eq. (|5.15|l . The data points come from the 
analysis of the Drell-Yan data from the E615 experiment [53l ]. 
The dashed line shows the recent reanalysis of the original 
data made in Ref. I70II. 



E615 Drell-Yan experiment [53]. The model results are 
represented with a band, which reflects the uncertainty 
in the determination of the scale Qo m Eq- (|5.15|) . The 
quality of this comparison is impressive, which shows that 
despite the rather embarrassingly low value of the scale 
Qo, the quark- model initial condition leads to fair phe- 
nomenology. The NLO evolution leads to small changes 
as compared to the LO results 30], in fact compatible 
with the experimental uncertainties. The dashed line in 
represents the recent reanalysis of the original E615 data 
made in Ref. [7(3]. We note that this result is also close 
to the band generated by our quark-model calculations. 

Moments of PDF's have been calculated on Euclidean 
lattices yielding (x) v = 0.3(1), (x 2 ) v = 0.10(5) and 
(x 3 ) v = 0.05(1) at Q = 2.4GeV. In Ref. [H (x) v = 
0.243(21) at Q = 2GeV. 

The non-singlet PDF in the pion was also evaluated 
on the transverse lattice [Til ] at the low renormalization 
scale Q ~ 0.5 GeV. In Fig[9]we show our PDF evolved to 
that scale (darker band) and to a lower scale of 0.35 GeV 
(lighter band). We take the liberty of moving the scale, 
as its determination on the lattice is not very precise. As 
we see, the agreement is qualitatively good if one con- 
siders the uncertainties of the data, especially when the 
lower scale is used. We also show the GRV98 parame- 
terization [72I ] (dashed line), which gives somewhat lower 
PDF (except low values of x) compared to the lattice 
data and our model results. 



, (°) 



-2C F 



n+l 



(n+l)(n + 2) 



k=l 



(5.17) 



with C F = (N 2 - 1)/(2N C ). The evolution equations 
can be solved via the inverse Mellin transformation. In 
Fig. [5] we confront the result for xq(x, Q) at the scale 
Q = 2 GeV with the data at the same scale from the 



D. Evolution of PDA 



PDA's have been intensely studied in the past in sev- 
eral contexts (see e.g. [ill, El [zl, [Zf3] and Ref. [zz| for a 
brief but comprehensive review). The PDA of the pion 
[3l| . which can be related to the isovector GPD through 
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xq(x) 




FIG. 9: (Color online) The quark-model prediction for the 
valence parton distribution function (PDF) of the pion for a 
single quark (either u or d for 7r + ) evolved to the scale Q = 
0.5 GeV (darker band) and Q = 0.35 GeV (lighter band). The 
width of the bands indicates the uncertainty in the initial scale 
The data come from the transverse lattice 
and correspond to the scale ~ 0.5 GeV. The 
parameterization at Q — 0.5 GeV. 



Qo, Eq. (15T5I 
calculations [7 
line shows the GRV [7 



the soft pion theorem, Eq. (|2.18|) . is 4>(x; Qo) = 1, which 
holds at the quark model scale Qo (|5.15|) 31]. The 
evolved PDA can be expressed in terms of the Gegen- 
bauer polynomials [zl, 



x,Q) = 6x(l-x)J2 C 3 J 2 {2x - l)a n {Q), (5.18) 

n=Q 



where the prime indicates summation over even values of 
n only. Our inial condition yields [3l| 



a n {Q) = 



2n + 3 



a(Q) 

3(n+l)(n+2) \a(Q ) 



7i°V(2/3o) 



,(5T9) 



where the anomalous dimension for the axial- vector ver- 
tex is the same as for the vector vertex, jn , given in 
Eq. (|5.17p . The evolved PDA is shown in Fig.[lj}l where 
it is compared to the E791 di-jet measurement [54| . The 
normalization of the di-jet data is used as a fit param- 
eter. Besides this normalization the result is parame- 
ter free. As we see the agreement with the E791 data 
is rather reasonable with a x 2 /DOF = 1.45. Nonethe- 
less, the asymptotic wave function generates a yet better 
X 2 /DOF = 0.45. Note that in our scheme such a extreme 
limit would correspond to taking Qo = Aqcd- 

The second Gegenbauer moment at the scale Q = 
2.4 GeV is ci2 = 0.12 to be compared with the value 
<i2 = 0.12(3) based on the analysis of the CLEO data of 
Ref. 80] where it was assumed that a n = for n > 4. 

Further, the leading-twist contribution to the pion 
transition form factor is, at the LO in the QCD evo- 




FIG. 10: (Color online) The quark-model prediction for 
the pion distribution amplitude (PDA) evolved to the scale 
Q — 2 GeV (band) and compared to the E791 di-jet mea- 
surement [Hit ] after proper normalization of the data. The 
width of the band indicates the uncertainty in the initial 
scale Qo, Eq. (|5.15p . We also show the the asymptotic PDA, 
4>(x,oo) = 6x(l — x) (dashed line). 



lution [78j, equal to 



Q 2 R 



7*- 



+ 7T7 (Q) 



2/ 



twist — 2 



6x(l — x) 



(5.20) 



The experimental value obtained in CLEO [8l| for the 
full form factor is Q 2 F 7 » i7r7 (Q)/(2/) = 0.83 ± 0.12 at 
Q 2 = (2.4GeV) 2 . Our value for the integral, 1.25 ± 0.10, 
overestimates the experimental result by about two stan- 
dard deviations, but one should bare in mind that higher 
twist contributions as well as NLO perturbative correc- 
tions have been ignored. 

The second ^-moment (£ = 2x — 1), defined as 



dx())(x,Q)(2x - l) 2 



(5.21) 



has been computed on Euclidean lattices yielding 
(£ 2 ) = 0.286(49) |3, (£ 2 ) = 0.269(39) H, (£ 2 ) = 
0.278(26) [HI from recent Euclidean lattice calculations 
at the scale Q = 1/a ~ 2.6 ± OTGeV, where a is the 
lattice spacing. Note that the asymptotic PDA would 
yield (£ 2 ) — 1/5 = 0.20. We get from the quark model 
(£ 2 ) = 0.244(4) for that scale. 

In Fig. [TT] we compare our model prediction for the 
PDA (band) to the transverse lattice data [7l[ at the 
scale Q = 0.5 GeV. A good agreement is observed. 

One of the most surprising aspects is that many of 
these results can be obtained from the integral rela- 
tion between PDF's and PDA's at a given scale, estab- 
lished in Ref. [3l|. The relation allows to predict <fi(x, Q) 
from V(x, Q) as parameterized, e.g., by the Durham 
group [62J. The method works for quark models, where 
V(x, Qo) = <j){x, Qo) = 1 at some scale Qo- 
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FIG. 11: (Color online) The quark-model prediction for pion 
distribution amplitude (PDA) evolved to the scale Q = 
0.5GeV (band) compared to the transverse lattice data [7l| . 
corresponding to the scale ~ 0.5 GeV. The width of the 
band corresponds to the uncertainty in the initial scale Qo, 
Eq. (I5T51 . 



E. End-point behavior 

The results at the quark-model scale exhibit disconti- 
nuity at x = 0, 1, as V(x, Qo) = 4>{ x i Qo) = 1- An impor- 
tant feature of the evolution is that it cures the end-point 
behavior of the PDF's and PDA's [H EE El • Usin g the 
Mcllin-momcnts formulation of the LO DGLAP evolu- 
tion it can be shown that if V(x, Qo) ~ c(l — x) p near 
x = 1 then 

V(x,Q)~c(l~x) p+8r ( Qo ' Q \ x->l, (5.22) 



V(x) 
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FIG. 12: (Color online) Comparison of the full result of the 
QCD evolution (solid lines) to the asymptotic formulas near 
the end-points (dashed lines). Top - valence PDF, bottom - 
PDA. The curves from top to bottom are for Q = 0.5, 2.4, 
and 10 GeV, respectively. The initial scale Qo is taken from 
Eq. (j5TT5]l . 



where we have introduced the short-hand notation 

KQo,Q) = ^log^ (5.23) 
2fj a{Q) 

and c and p are some constants. Specifically, in our 
quark-model case c = 1 and p = 0. The prefactor can 
also be obtained, yielding the more accurate formula (8o| 



V(x,Q) 



e 2(3-4 7 )r(Q ,Q) 

r(l + 8r(Q ,Q)) 



(1-x) 



Sr(Q ,Q) 



(5.24) 



where x — ► 1, 7 denotes the Euler-Mascheroni constant, 
and r is the Euler gamma function. At Q > 1 GeV the 
exponent of 1 — x is a function weakly dependent on Q. 
The explicit forms for several values of Q and Qo from 
Eq. (f5TT5|) are 



V(x,0.5 GeV) 
V(x,2A GeV) 
V(x, 10 GeV) 



1.23(1 
1.13(1 
1.00(1 



0.0.3 



x) 
x) 1A5 



(5.25) 



Similarly, for the PDA evolved with the help of the 
Gegenbauer polynomials (see Appendix [E] for the deriva- 
tion), one can show that 5 



<j){x,Q) 



r(-4r(Q ,Q)) 2(3 -4 7 )r(Qo,Q) 

r(4r(Q ,Q)) 



(5.26) 



and a symmetric expression when x — > 1 with x replaced 
by 1 — x. With several explicit values for Q we have 



„0.26 



(x,0.5 GeV) - 1.51 x u 
(x,2A GeV) - 2.95 x 59 , 
i(ar, 10 GeV) - 4.65 x 73 , 



(5.27) 



all for Qo from Eq. (|5.15[) . These asymptotic forms are 
compared to the full result of the QCD evolution in the 
bottom panel of Fig. [T2j We note that the range of valid- 
ity of the approximation (|5.26|) shrinks closer and closer 
to the end point as Q is increased. This must be so, 



These are compared to the full result in the top panel of 
Fig. [TU 



5 The corresponding formula in Ref. [3lll has a mistake. 
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as at Q — > oo the asymptotic form 62 (1 — x) sets in in 
the whole range of x, while the power r(Q , Q) increases 
indefinitely with Q. 

It should be noted that the behavior of 
Eqs. (|5.22l5.26p . exhibiting the desired continuity 
of the functions at the end-points, is achieved already 
at values of Q infinitesimally larger than Qq. Thus the 
QCD evolution heals the end-point problem immedi- 
ately, at any Q > Qq. Such a phenomenon is linked to 
the nonuniform convergence of the Mellin or Gegenbauer 
functional series near the end-points for the PDF and 
PDA, respectively. 



F. Evolution of diagonal GPD in the 
impact-parameter space 

The impact-parameter dependence quoted in 

Eq. (|3.14l) at the quark model point [32| not only 
satisfies the model independent relation (|2.15[) but 
after proper smearing over plaquettes and DGLAP 
evolution qualitatively reproduces both the Bjorken x 
and impact parameter dependence when compared to 
transverse lattice results [H, |87} at the rather low scale 
Q ~ 0.5GeV. This is a remarkable finding, since the 
transverse lattice at such low scales should incorporate 
non-perturbative evolution effects if they happened to 
be important. Details can be found in Ref. [32j. 



G. Discussion 

To sumarize this Section, the low value of the renor- 
malization scale Qq deduced from the LO perturbative 
evolution of the momentum fraction, complies surpris- 
ingly well with a wealth of fragmentary information for 
the non-singlet partonic distributions both on the exper- 
imental side as well as compared to Euclidean lattices at 
Q ~ 2GeV and transverse lattices at Q ~ 0.5GeV. This 
provides some confidence on applying a similar strategy 
to the evolution of non-singlet GPD's as we do in the next 
section. Of course, it would be of great help to have Eu- 
clidean lattices at small renormalization scales, such that 
some of the non-perturbative evolution could be explic- 
itly seen. Unfortunately, the transition from the intrin- 
sically non-perturbative lattice regularization to the per- 
turbative MS regularization scheme requires high scales, 
so such a calculation seems hardly viable. Transverse lat- 
tices do not suffer from this drawback, as these are non- 
perturbative calculations at low scales [ll|. We are in 
qualitative agreement also with these lattice calculations, 
which probe the evolution in a region where it might po- 
tentially be highly non-perturbative. Our analysis agrees 
within uncertainties with a picture where the main non- 
perturbative feature of the valence quark contribution is 
provided by the initial condition. In any case, as shown 
in the LO and NLO analysis of Ref. [30] , the sea-quark 
and gluon PDF's from chiral quark models are less prop- 



erly reproduced. This might be improved if some non- 
singlet either sea or gluonic model contributions could 
be provided at the model scale, Qq. Despite the efforts 
all over the years the problem of determining the non- 
perturbative gluon content in a hadron at low scales has 
remained unresolved. These provisos should be taken in 
mind when evolving the singlet GPD's in our scenario. 



VI. THE QCD EVOLUTION OF GPD'S 

The explicit form of the LO QCD evolution e quat ions 
for the GPD's can be found in [|| M, M, M, M, M [HI . 
In this paper we solve them with the numerical method 
developed in [§2], based on the Chebyshev polynomial 
expansion. 

As extensively discussed in the previous section, per- 
turbative QCD brings in the issue of evolution equa- 
tions for the GPD's. Similarly to the more familiar 
case of the PDF, the QCD interactions of massless par- 
tons lead to collinear divergences which are factored out 
and absorbed into the GPD's. As a result, the GPD's 
become dependent on a factorization (renormalization) 
scale /i, usually identified with the hard scale, fi = Q. 
Thus, in general, the GPD's are functions of four vari- 
ables, H = H(X, £, t, Q 2 ), with the kinematic constraints 
|A| < 1 and < £ < 1. The renormalization group 
equations which govern the dependence of GPD's on Q 2 
are described in detail e.g. in Ref. [92]. The form of 
these equations depends on the asymmetry parameter 
£, which defines two regions: the Efremov-Radyushkin- 
Brodsky-Lepage (ERBL) region for |A| < £, and the 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) re- 
gion for |X| > £. 

An important feature of the GPD evolution, which 
makes it more complicated than in the case of PDF or 
PDA, is that the evolution equations in the ERBL region 
depend on the values of GPD's in the DGLAP region. 
The converse is not true, the evolution in the DGLAP 
region is not influenced by the the ERBL region. 

Asymptotically, for Q 2 — > oo, the GPD's tend to the 
asymptotic forms which are concentrated in the ERBL 
region only. In particular, for |X| < £ we have 



H 



i=i 



(6.1) 



H 



1=0 



15 



N 



f 



X ( A 2 ' 



xh 9 = a-e 



15 



4C F 



Lf>;4fy + \/ v 1 e ' inn 



while the GPD's vanish for |A| > £ . The proportionality 
constants reflect the normalization of the GPD's at the 
initial scale Qo, as the following charge- and momentum- 
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FIG. 13: (Color online) Results of the LO QCD evolution from the SQM initial condition for several values of t and £ = 1/3. 
Solid - initial condition at the quark-model scale, dashed - evolution to Q 2 = (4GeV) 2 , dotted - asymptotic form, Q 2 — > oo. 
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FIG. 14: (Color online) Same as bottom left panel of Fig. [13] 
for the NJL model. 



conservation sum rules are preserved by the evolution 
J dXH I = 1 (X,£,t,Q 2 ) = 2F v (t), (6.2) 
J dX (XH I=a (X, e, t, Q 2 ) + XH g (X, e, t, Q 2 )) 

in accordance to Eq. (|2.9I2.10[) . 

The results of the LO evolution from the SQM ini- 
tial condition at the scale Qo up to Q = 4 GeV and 
£ = 1/3 are shown in Figs. [T31 In each set of the four 
plots, corresponding to a different value of t, we show the 
quark non-singlet H I=1 (top) and singlet (middle) H I=a 
distributions together with the gluon H g (bottom), con- 
ventionally multiplied by X. We have chosen the sample 
value £ = 1/3, since the results are qualitatively simi- 
lar for other values of £. The solid lines show the initial 
condition at the quark- model scale Qo of Eq. (|5.15[) , the 
dashed lines show the result of the LO QCD evolution 
to the scale Q = 4 GeV, and the dotted lines show the 
asymptotic forms at Q — > oo given in Eq. (|6.ip . As the 



FIG. 15: (Color online) Same as Fig. [13] for subsequent evo- 
lution scales: Q 2 = 0.1, 1, 10, 10 2 , . . . , 10 s GeV 2 . Higher Q 2 
gives higher magnitude of the curves in the ERBL region. 



value of —t is increased, the magnitudes of the curves be- 
comes lower, conforming to the sum rules (|2 .912 . 1 0|) . We 
note that the evolution smooths out the original distribu- 
tions, in particular, the discontinuities at the end-points, 
X = ±1, and at the ERBL-DGLAP matching points 
X = ±£ disappear for the isosinglet GPD. 

The results for the NJL model are very similar to the 
case of SQM. In Fig. [M] we show them for t = -1 GeV 2 
and £ = 1/3. This similarity between the models is a 
sheer reflection of the numerical similarity in the initial 
condition, cf. Fig. QJ for t = -1 GeV 2 and Fig. [HI 

In Fig. [T5] it is shown how slow the evolution is in 
reaching the asymptotic forms of the GPDF's. The evo- 
lution is fastest at low values of Q, where the coupling 
constant is large, and it immediately pulls down the end- 
point values to zero. Then, the strength gradually drifts 
from the DGLAP regions to the ERBL region. Yet, the 
approach to the asymptotic form is very slow, with the 
tails in the DGLAP region present. The highest Q 2 dis- 
played in the figure is 10 8 GeV 2 and the asymptotic form 
is reached at "cosmologically" large values of Q, which 
are never achieved experimentally. Thus, the only way to 
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approach the asymptotics would be to start from initial 
conditions which are already close to it. We also observe 
that a larger —t more strength of the quark GPDs resides 
in the DGLAP region. This feature reflects the shape of 
the initial condition, which inhibits the strength in the 
ERBL region. Evolution up to Q = 2 GeV retains this 
behavior, which gradually disappears as Q — > oo where 
all the strength settles in the ERBL region. 

We note that the desired vanishing of the GPD's at 
the end-points X = ±1 is achieved due to the QCD 
evolution, similarly to the results presented in Sec. IV El 
Also, evolution leads to continuity at the DGLAP-ERBL 
boundary, X = ±£. These features are achieved at scales 
Q infinitcsimally above Qq. 



VII. CONCLUSIONS 

We summarize our main points. In the present paper 
we have dealt with the determination of the leading twist 
GPD's of the pion in field theoretic chiral quark models. 
We have done so with the help of an efficient method 
using the a-representation of the quark propagators and 
an extensive use of double distributions. Our calculation 
incorporates the necessary D-terms required by polyno- 
miality and dimensional analysis. In the chiral limit, we 
have been able to determine explicit analytic formulas 
for the pion GPD's. All a priori properties which ought 
to be satisfied on general principles, namely polynomial- 
ity, positivity, proper support, soft-pion theorems, sum 
rules and normalization arc indeed fulfilled explicitly by 
our model calculation. Although one might superficially 
think that these properties should be trivially satisfied, 
the fact that one deals with regularization or momen- 
tum dependence makes the fulfillment of those properties 
less obvious, and in fact many calculations violate some 
property. A key ingredient in our approach has been 
a scrupulous treatment of regularization in conjunction 
with electromagnetic and chiral Ward-Takahashi identi- 
ties. Our results for the pionic GPD's in the NJL model 
agree with Ref. 

In the two chiral quark models considered, NJL and 
SQM, we have found results looking alike since the mod- 
els are mainly fixed by the pion charge form factor which 
in both cases looks very similar. In addition, we have 
determined the pion gravitational form factor entering 
the momentum-conservation sum rule. The outcoming 
GPD's arc not t-factorizable, an assumption which is be- 
ing extensively used in phenomenological approaches, not 
based on consistent dynamical calculations. 

However, with all those desirable properties fulfilled, 
one must undertake the ERBL-DGLAP QCD evolution 
in order to relate the model results to experimental 
high-energy data where higher-twist contributions to the 
GPD's can be disregarded. Likewise, a comparison to 
lattice results of the twist-2 GPD's requires specification 
of a running scale. This aspect of the calculation is most 
frequently ignored in dynamical model calculations, and 



particularly in chiral quark models. A practical compari- 
son to either experiment or lattice can be be achieved by 
matching the momentum fraction of the QCD evolved 
quark model to the experimental or lattice-extracted re- 
sult. In practice, the LO perturbative evolution is used 
with the result that the low-energy quark-model scale 
is very low. Nevertheless, once this is fixed the GPD's 
are uniquely determined. We have confronted our pre- 
dictions with all available information extracted either 
from experiment or lattice both transverse or Euclidean. 
The experimental data include the Fermilab E615 and 
E791 measurements of PDF and PDA of the pion, re- 
spectively, and CLEO measuraments on the pion transi- 
tion form factor. The reasonable overall agreement to all 
these data corresponding to quite different kinematical 
situations should be stressed. This fact provides some 
confidence on our predictions of the non-singlet lead- 
ing twist GPD's. This also applies to the pion-photon 
Transition Distribution Amplitude, determined recently 
in quark models [35| . 
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APPENDIX A: THE a-REPRESENTATION 
EVALUATION OF THE TWO- AND 
THREE-POINT FUNCTIONS 

In studies based on quark models the simplest way to 
obtain the GPD's and the double distributions is through 
the use of the a-representation for the scalar propagators. 
The advantage of this representation over other popular 
Feynman parameterizations of one-loop functions relies 
in the fact that the S(k ■ n — x) function constraining the 
loop integration is also naturally written in terms of an 
integral of an exponential. This allows for maintaining 
the explicit Lorentz covariance throughout the calcula- 
tion. In our scheme one does not have to start with the 
somewhat cumbersome moments in k ■ n, and then "in- 
vert" the result, as is frequently done. The method used 
in this paper leads to well-defined and very simple algebra 
and reproduces the double distributions in chiral quark 
models from the literature. Also, in our approach the 
otherwise subtle effect of the emergence of the D-terms 
follows in a clear way just from the Feynman diagrams. 

Below we derive basic integrals appearing later on in 
the evaluation of the GPD's. Calculations of this Ap- 
pendix are made in the Euclidean space. We denote the 
Euclidean scalar propagators of particles of mass lo as 

Sk = TT = WT12 = / dae~^ 2 ^ 2 \ (Al) 
D k k 1 + uj z J Q 

where the RHS displays the a-representation. 
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1. Two-point functions This two-point function with the constrained k ■ n inte- 

gration can be written as 

Let us first consider the function I(x, I ■ n, 0, I 2 ) corre- 
sponding to the definition (|2.26|) with the choice I' = 0. 



I[X,K, 0, I ) = / /O-m "^'"'^^'^-^ - ^ 



P J (2^) 4 
AN c w 2 f d 4 k' 

AN r w 2 f d 4 k' 



d 4 k 

WT 4 



(2tt)< 



OO pOC 

da / db 
i' Jo 

OO pOO 



dX 

b j2 



OO poo 



e iX(k-n-x) da dbe -a{k'+^)-H(k-lY+^) 

Jo Jo 



(a + b)e~ {a+h){k ' 2+u2) ~^ l2+tX ' [hK ~ (a+V)x) 



da / dbS[bK- (a + b)x}(a + b)e- {a+b){k ' 2+uj2) -^ l \ 



(A2) 



where for brevity k — n ■ I, the shifted integration mo- 
mentum is k' = k+ 7 ± T l + ±\'n, and A' = X/(a + b). The 
S [bn — (a + b)x] function gives the constraint 



b' 



(A3) 



Since the integration variables are positive, a, b > 0, 
it follows immediately that x <G [0,n] for n > and 
x G [— k, 0) for k < 0. This provides the proper sup- 



port for I(x, k, 0, 1 ), which can be written generally as 
9[x(k — x)]. 

One can decompose the k' integration into two parts, 



dk m dk' 3 = ndK", dk' l dk rz = ndu, 



(A4) 



with K 2 = (k'°) 2 + (fc' 3 ) 2 and the "transverse" momen- 
tum u = (k' 1 ) 2 + (k' 2 ) 2 . Then 



I{X,K,0,1 2 ) 



N c uj 2 9[x(k - x)] 

4tt 2 / 2 
N c w 2 0[x{k - x)] 
4tt 2 / 2 |k| 



OO pOO poo 



du da dbS[bn- {a + b)x\e- [a+b)[u+UJ } ~*^ 
i! Jo Jo 

/ -f,'[«+«. 2 +f(i-f)] = N c u 2 6[x{k - x 



(A5) 



OC P OO 

du I db' e 
o Jo 



du 



Aw 2 f 2 \ K \ J u + w 2 + z(l-z)p> 



where b' — bn/x. The integral over u is logarithmically 
divergent, hence needs regularization, as expected. 

The general function I(x,k,k',(1 — I') 2 ), where 
k' = n ■ I', involves no extra work, as it can be obtained 
from the I' = case with the replacement 

k — ► k — I , I — > I — I , x — > x — k , k — > k — k . (A6) 

This yields 



I{x,k,k',(1-1') 2 ) = 



du 



N c oj 2 6[(x- k')(k-x)} 

4lT 2 f 2 \K - K'\ 



(A7) 



An important consequence of Lorentz invariance is 
polynomiality [l[ . We verify it by introducing the variable 
v = (x — k')/(k — K r ), when I(x, K, n' , (I — I') 2 ) becomes 
a function of v devided by \k— k'\. We obtain (assuming 
for definiteness k > n') 



dxl(x, k, k', (I — I') )x n — 



dvf(v)[K' + v(k- n')] n , (A8) 



which results in a polynomial in n and n' of the order at 
most n. The first few moments have the explicit form 
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i V ' ' ' ' 4tt2/ 2 J ' 
dxl[x, k, k ,t) x = — — j-Tq- / aw- 



-i v ' ' ' ' ^ 2 P 7o 

axlyx, k, k ,t) x =— — 5-T5- / dux 
-i 4tH/2 y 

(4A + 1) (« - K ') 2 - ViA+Tlog ( ^jggl ) ((2^ + 1)« 2 - 4Aai'/c + (2A + 1) (nf) 

AA+ 1 



(A9) 



where A = (it + w 2 )/r and r = (I - I') 2 . 

In the literature the D-term is by definition the two- 
point function in the i-channel . It originates from the 
diagram with the contact pion-quark term as well as from 
the reduced three-point diagram, where by "reduction" 
one means the replacement of k 2 and k-l pieces appearing 
in the numerator from the trace factor, in terms of the 
inverse scalar propagators. The two-point functions in 
the s-channel (resulting from the reduction of the three- 
point function) are traditionally treated as singular parts 
of the double distributions. 



2. Three-point functions 



For the three-point functions we proceed analogously, 
now with three scalar propagators. We need to take into 
account the kinematics of the direct and crossed diagrams 
of Fig. [TJ We first analyze in detail the the three-point 
function resulting from the direct diagram (a), since the 
case of the crossed diagram is obtained via a simple kine- 
matic transformation. We have 



J(x,q-n,p-n,q 2 ,p 2 ,p-q) = J -^-^5(k ■ n ~ x)S k S k+q S bk - p (A10) 

= AN c u 2 f ( [ -e 1 ^^' 1 ' [la [°db [ d ce --(fcW)-K(fc+9)W)-c(( fc -p) W). 



Shifting the integration variable, k' — k + (f3q — "fp — iXn/2) /{a + + 7), and carrying over the d 4 k' integration yields 

j = r da r db r dc 1 ,L ^"-^" ) e -(^-^^^-^p, (AU) 

4^ 2 / 2 7 Jo Jo (a + b + c) 2 V a + b + c J 



Next, we change the variables into 

b 

s = a + o + c, y = - , 



(A12) 



Note that since a,b,c> 0, we get < y, z < 1 and also 
y + z < 1. The substitution and integration over s yields 



J 



N c lo 2 
4tt 2 / 2 

[ dy 



dz- 



(A13) 

9(1 — y — z)5 (x — zp ■ n + yq ■ n) 
'■ + y(l — y)q 2 + z(l — z)p 2 + 2yzp ■ q 



We note that polynomiality is obvious from this form, 
as multiplication by the power x n is equivalent to the 
multiplication by the factor (zp ■ n — yq ■ n) n . In the 
chiral limit of m n — the first few moments are relatively 
simple: 
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arctan 



-2 i Vt 



sjAw 2 - t 



arctan 



Jo 



<l.r -fx 2 = N n c ^.„ [t(C(C + 3) - 7) + 



arctan 



47T 2 / 2 * 2 
V4u> 2 - t 



{^/ty/Aw 2 - t(6 - C(C + 2)) + 2 (* - 2w 2 (C - 1)) arctan ( ~/=f ) 



(A14) 



(A) ] 2 - V ^ (C ' 2) arCtan (A) + ^ - ^ ' 



We can rewrite Eq. (|A13[) as 

N c lu 2 9(1 - y - z) 



J(x) 



4tt 2 / 2 uj 2 - y(l 
i ,i 

dy / dz5 (x - 



- y — z)t — z(\ — z)m 2 ' 

(A15) 



where we have used the kinematics (12.11) . The curly T 
denotes the double distribution. 
Let us denote 



V 



2/(1 - V- z)t - z(l - z)m\ 



(A16) 



For the GPD of the pion, due to the crossing symmetry, 
one may assume < C < 1- Next, we perform the z 
integration, which sets 



(A17) 



The distributions in Eq. (| A13|) give the following limits 
for the y integration: 



J = 



N c uj 2 
An 2 ! 2 



1[x(C-x)] +9[(x-C)(l-x)] 
Jo Jo 



1-x ' 

l-C 



dy 
) V 1 

(A18) 



with the first term having the support x £ [0, £], and 
the second x e [C, 1] - The function F(x) is continu- 
ous, but the derivative dF(x)/dx is discontinuous at the 
points x = 0, C 1. The double distribution is Miinchen- 
symmetric [95[, i-e. T(z,y) = T(z,l — y — z). This 
feature is related to the crossing symmetry, holding for 
identical particles. 

The result for the crossed diagram (see Fig. [T] is ob- 
tained from the above result for the direct diagram with 
the replacement p — > — p — q. Replacing correspondingly 
x — > C — x and performing the Miinchen transforma- 
tion [9a 



y 



i 



y 



(A19) 



we find that T> is invariant under this joint transforma- 
tion. The function 5(x — z — y£) is also invariant under 



these combined two trasformations. Finally, the support 
is invariant, since 

9(l-y-z)9[y(l-y)]9[z(l-z)}^ 
9(y)9[(l-y-z)(y + z)]9[z(l-z)\ 
= 9(1- y- z)9[y(l - y]9[z(l - *)], (A20) 

where the equality in the above formula is an algebraic 
identity. Therefore the crossed diagram is related to the 
direct diagram as follows: 



crossed 



Or,C) = J dircct (C-^C) 



(A21) 



The support of the crossed diagram reflects the support 
of the direct diagram, i.e., x G [—! + £,£]. 



APPENDIX B: THE TWO- AND THREE-POINT 
FUNCTIONS IN THE SPECTRAL QUARK 
MODEL 

According to the general rule, in SQM one append the 
formulas with the spectral integration J c dojco 2 p(oj) . The 
results below are for the meson-dominance model. 



1. The two-point function 

We assume k' < K. The spectral integration yields 



Isqm(x,k,k ,1 )= / dww p(u)I(x,k,k ,1 
Jc 

9[(x - k')(k - x)} 

3/2 • 



where we have used the relation 

My = 2An 2 f/N c . 
The integration over x yields the form factor 

dx!sQM(x, K, k', I 2 ) 



M 2 



M 2 + I 2 ■ 



(Bl) 



(B2) 



(B3) 
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In agreement with polynomiality, this form factor is in- 
dependent of the value of k or n' . Note that the vector- 
meson dominance is readily obtained. Similarly, 



dxL 



SQM 



M v k + k' 
M 2 +l 2 ^2~' 



(B4) 



2. The three-point function 

For simplicity in this Appendix we work in the chiral 
limit. In this case the double distribution becomes 



Nm 2 



(1 -!/-*) 



4tt 2 / 2 - y(l - y - z)V 



(B5) 



and the subsequent spectral integration yields (again we 
only multiply by u> 2 and leave out other factors) 



FsQM{z,y;t) = J dwp(u))F u {z,y;t) 

0(1- y-z) 



c 

Nr. 



47T 2 / 2 / _ 4y(l-y-z)t \ 5 / 2 

1 Ml 



Let us introduce the short-hand notation 



(B6) 



Xi 



Xx 



2{x - 1) [3(C - 1)M 2 + t(x - l) 2 ] 

[(C-l)A4 + ^-l) 2 ] 2 ' 
(x(C - 2) + C) (3M 2 (C - 1)C 2 + * ((C 2 + 8C - 8) x 2 + 2(4 - 5C)C* + C 2 )) I 



((C 



1)M 2 



t(x-l) 2 ) 2 (c 2 



4tx(x-C) 



3/2 



+ oX2- 



(B7) 



Then, for the case £ > 0, 

J SQ u(x,C;t) = (0[x(C-x)] X i + 0[(i-^)(x-C)]x2). 

(B8) 

The function satisfies Jsqm(0, Ci*) = ^SQm(LC^) = 0- 
The value at the matching point x = ( is 



which is a triangle of area 3 /M v = N c / '(8tt 2 f 2 ). 
For the case C, = we have 

(3M 2 - til - x) 2 ) (1 - x) , 

JsQMfe0;!) " -.)»)■ ' (B13) 



•^sqm(Ci C; *) 



2 (3M 2 + t(C - 1)) 
(M 2 +t(C- I)) 2 



(B9) 



APPENDIX C: THE TWO- AND THREE-POINT 
FUNCTIONS IN THE NJL MODEL 



The integration over x produces a ^-independent (as re- 
quired by polynomiality) form factor, 



dxJ SQM {x,C,t) = M 2 _ I 



(BIO) 



Operationally, the calculation in the NJL model with 
the regularization (|4.1[) amounts to taking the generic 



expressions (|A7j) and (|A18|) . replacing J 1 — > M 2 + A 2 in 
the denominators, carry out the integrations, and finally 
applying (14. ip . We work in the chiral limit. 



Similarly, 



/' 

Jo 



dxJsQM.(x,C;t)x 



l0g i 1 ~ 



My ~ t 



t 



(Bll) 



For the special case of t = Eq. (|B8|) reduces to the 
very simple expression 



+ 9[(l-x)(x-Q] 



x- 1 



(B12) 



1. The two-point function 

Through the use of Eq. (|4.2p we arrive immediately at 
the formula (for k> n') 



/njl(x,k,k',0) 



0[(x - k')(k - x)] 



(CI) 



2. The three-point function 



We find 
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Jnjl(x, C;"V = 0; A) = 
log 



N r M 



8tt 2 / 2 

~{x-l) + ^ -t(x-l) 2 -A(Q-l){M 2 +\ 2 ) 



2 log 



\/^t( a: -l) + ^/-t( a: -l) 2 -4(C-l)(Af 2 +A 2 T 
^/-t(x-l) 2 -4(C-l)(M 2 +A 2 )-\/=t(ii;-l) 



0((l-aO(s-O) 



v / -t(x-l) 2 -4(C-l)(M 2 +A 2 )- v ^t(i:-l) 



log 



t^/-t(a: - l) 2 - 4(C - 1) (M 2 + A 2 ) 

/ y / -t(x-l) 2 -4(C-l)(Af 2 + A 2 )C+v^t( a: (C-2)+C) 

\ Cv / - t ( 3; - 1 ) 2 - 4 (C-i)(Af 2 +^ 2 )-%^t(^(C-2)+C) 



(C2) 



0(a;(C - x)) 



ty/-t(x - l) 2 - 4(C - 1) (M 2 + A 2 ) 



and 

^njl (fj C; ™7r = 0) = Jnjl(x, C; w w = 0; A)| rcg . (C3) 

Polynomiality follows from the fact that the expressions 
are derived from the double distributions (|A15j) and 
the distributive nature of the regularization prescription 
(|4.1[) . By distributive we mean that it is a sum over quark 
masses or the integral over u of the formal expressions 
for / and J. 



APPENDIX D: THE GRAVITATIONAL FORM 
FACTORS 

The gravitational form factors of the pion (5f| are 
defined through the matrix clement of the energy- 
momentum tensor, 



(7^ + 9)10^(0) |tt») 

S ab 



(Dl) 



^ OTY - (q^nOiiq 2 ) + (2p + qY(2p + q) v 6 2 (q 2 )) . 

They satisfy the low energy theorem 9i(0) — 62(0) = 
0(m 2 ) [Hll]. The leading- N c quark- model evaluation 
amounts to computing the diagrams of Figs. [Hand [2] with 
the pion gravitational vertex 



9^(k + q,k) = - ((2k + q)»>f + (2fc + q) V) 



(D2) 



The results of the calculation in SQM is Eq. (|3.5p . Equa- 
tion (|2.10[) follows from considering the matrix clement 



oin^ v n v . Then 

(n b (p + q) \n^(0)n v | tt») = (D3) 
S ab \[C 2 9i(q 2 ) + (2-0 2 d2(q 2 )]. 

The vertex becomes 

n„e" v (k + q, k)n v = (x - C/2) 7 • n. (D4) 

We notice it is the same vertex as in the evaluation of 
the GPD's multiplied by (x — C/2). Upon passing to the 
symmetric notation Eg. 12.101 follows. 



APPENDIX E: END-POINT ANALYSIS FOR THE 

PDA 

Here we derive the formulas used in the main text 
of Sect. IV El for the PDA and correct a mistake in ex- 
pressions of our previous work [2?], Right at the 
end-points, x — 0,1, the series (|5.18j) diverges since 
C^ 2 (±1) = |(2fe + l)(2fc + 2), meaning a non-uniform 
convergence as x — > or x — * 1 as well as the large-n 
dominance of the end-point behavior. In this limit we 
have IH 



a{Q) 
a{Qo) 



7™/(2/3o) 



n -8r(Q ,Q) e 2(3-4 7 )r(Q ,Q)^ 



(El) 



Only even-n terms contribute in Eq. (|5.18j) . hence we 
impose this condition for integer n and then extend n to 
real values. The summation in Eq. (|5.18|) is then replaced 
with an integral, 

1 r°° 3 

4>(x,Q) ~ -8a;e (3/4 - 7)Q / dn C„ 2 (2x - l)n- a_1 , 

2 Jo 

(E2) 

where the factor of 1/2 comes from the summation over 
even n only, and 



a = 8r(Q ,Q). 



(E3) 



The Gegenbauer polynomials C„ (£) in the variable 
£ = 2x — 1 satisfy the differential equation [96| 

(1 - e)y"(0 - 4£y'(0 + n(n + 3)y(0 = , (E4) 
which upon the substitution 

«(0 



y(0 



i-e 



(E5) 



transforms into a Schrodingcr-like equation at zero en- 
ergy, 



u"(0 



n 2 + 3n + 2 



u(0 = 0. 



(E6) 
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Here the interval of interest, — 1 < £ < 1, corresponds to 
the classically allowed region and the potential is attrac- 
tive. In the large-n limit the solution oscillates rapidly 
and a semiclassical WKB approximation [97| might be 
used. Here in order to analyze the limit x — * we con- 
sider the differential equation, Eq. (|E6j) . which in the 
limit x — > and n 3> 1 transforms into a zero-energy 
Coulomb-like problem. Its solution can be generally writ- 
ten in terms of the Bessel functions, 



Crl (2S - 1) 



CxJi{2n^fx) + c 2 Y 1 (2n^/x) 



(E7) 



where actually C2 = 0, since the Gegenbauer polyno- 
mials are regular at the end-point x — 0. The unde- 
termined constant c\ may be obtained by matching the 
small argument expansion J\{z) = z/2 + . . . to the value 

Cn^ 2 {{) ~ \, yielding for n ^> 1 and x — > the formula 



Ci{2x-l) 



2^fx~ 



Ji(2n\/x) . 



(E8) 



As we see in the \ow-x and large-n limit there is a scaling 
behavior of the Gegenbauer polynomials, 



C£(x) ~ nF{ny/x)/y/x . 



(E9) 



From Eq. (|E2p it is now clear that this translates into the 
low- a: scaling behavior 



poo 

(x,Q) ~ x a/2 e {3/4 ^ )a dtF{t)t- a . (ElO) 
Jo 



Evaluation of the integral and collecting the factors yields 
Eq. {PUD - 
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